2-d Self- Avoiding Walks on a Cylinder 
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We present simulations of self-avoiding random walks on 2-d lattices with the topology of an 
infinitely long cylinder, in the limit where the cylinder circumference L is much smaller than the 
Flory radius. We study in particular the L-dependence of the size h parallel to the cylinder axis, the 
connectivity constant n, the variance of the winding number around the cylinder, and the density of 
parallel contacts. While /j,(L) and (W 2 (L, h)) scale as as expected (in particular, {W 2 (L, h)) ~ h/L), 
the number of parallel contacts decays as h/L 192 , in striking contrast to recent predictions. These 
findings strongly speak against recent speculations that the critical exponent 7 of SAW's might be 
nonuniversal. Finally, we find that the amplitude for {W 2 ) does not agree with naive expectations 
from conformal invariance. 



PACS numbers: 05.70.Jk, 61.25.Hq 

Although self-avoiding random walks (SAW's) are 
among the best studied critical phenomena, there have 
been recently speculations which, if true, would have 
very far reaching and surprising consequences Q. One 
such consequence would be that the critical exponent 
7 of interacting SAW's on the square lattice would be 
temperature dependent if the interaction between neigh- 
boring bonds (i.e. steps of the walk on opposite edges 
of a plaquette) depends on the relative orientation of 
the steps. Another consequence would be that 7 for 
a-thermal SAW's on the Manhattan lattice is different 
from 7 on other 2-d lattices (^|. Numerical evidence for 
the former prediction is lacking ||-|]], but it has been 
argued Q that this is not conclusive since much longer 
SAW's (10 9 steps) would be needed to refute the predic- 
tion unambiguously. As concerns the Manhattan lattice, 
the numerical evidence is unclear since non-universality 
was found in but universality was found in |Io| . 

The main question at stake here is the density of par- 
allel contacts (i.e. the number of pairs of parallel steps 
on opposite sites of one plaquette) in a very long SAW. 
It is easily seen that parallel contacts are forbidden for 
2-d self-avoiding closed loops. Therefore, it seems plausi- 
ble that parallel contacts result only from spirals formed 
by the ends of the walks. This was indeed suggested by 
simulations in Q , but again it seems difficult to perform 
significant measurements on planar 2-d SAW's. 

As pointed out in the natural geometry to study 
this problem is obtained by mapping the plane onto the 
surface of a cylinder by the conformal map 



z = x + iy—>w = (L/2n) lnz. 



(1) 



Except near the ends of the walk, parallel contacts can 
only occur when the walk wraps around the cylinder. In 
the limit of very long chains (N u ~S> L, where N is the 
number of steps) a typical SAW has to grow either paral- 
lel or antiparallel to the cylinder axis. Finite size scaling 
predicts that its longitudinal size h (defined here as the 



average end-to-end distance, projected onto the direction 
parallel to the cylinder axis) scales as 



N/L 



3/4 



(2) 



The scaling of the winding number variance follows from 
the facts that it must be linear in h, and that it is di- 
mensionless, 



(W 2 ) « A h/L 



(3) 



where A should be a universal amplitude. In the fol- 
lowing, winding numbers will be measured in terms of 
wrappings around the cylinder, with W = 1 correspond- 
ing one full turn. Finally, the connectivity constant (the 
critical monomer fugacity) should show the usual finite 
size behavior 



(4) 



Here, the amplitude b is not universal. A universal am- 
plitude can be obtained for the free energy F = — In Z w 
N ln/i, where Z is the partition sum. Using eq.(||) we 
obtain 



F{L,N)=F oa {N)-Bh/L, 



(5) 



up to terms logarithmic in N . 

As discussed in Q , the predictions of [Q can be under- 
stood intuitively from the assumption that the density of 
parallel contacts follows essentially the winding number: 
If the walk wraps once around the cylinder, there should 
be O(l) parallel contacts. Thus the prediction for the 
average number of parallel contacts is 



C h/L 



(6) 



with C being yet another universal amplitude. 

Finally, we can rewrite eq.(|J) in terms of the Flory 
radius R of SAW's in planar geometry. More precisely, 



1 



we define R as the rms. 
SAW of iV steps, R~ N u . 



end-to-end distance of a free 
Then we obtain 



*=d(* 

L \L 



l/u 



(7) 



with D being also universal. Using this, we can eliminate 
h/L from eqs.(|),(]|), and (§) to obtain ((W 2 ), F^N) - 
F(L,N), n||) = (A',B',C) (R/L) 1 ^, with universal 
amplitudes A' = AD, B' = BD, and C = CD. 

The main objective of the present paper was to test 
eqs.(Q-^[) by means of Monte Carlo simulations. Since we 
want to simulate lattices with large perimeters L, and 
with N v » L (we went up to L = 128 and N = 60, 000), 
we need a fast algorithm. For SAW's in planar geometry 
the fastest known algorithm is the pivot algorithm (at 
least if one is not interested in a). But it is easy to see 
that the pivot algorithm doesn't work well on the cylin- 
der (for sufficiently small L one can prove that it isn't 
even ergodic). We therefore used the pruned-enriched 
Rosenbluth method (PERM) |ll|] with markovian antici- 
pation flj,[l3| • In PERM one starts off by growing chains 
according to the well known Rosenbluth method |TJ| , but 
when their weights become too large resp. too small one 
interferes by cloning resp. pruning. Weights have to 
be used since Rosenbluth sampling is biased, and the 
weights are needed to compensate the bias. 

In fc-step markovian anticipation one uses an addi- 
tional bias based on the statistics of sequences of k + 1 
successive steps. Let us number the 2d directions on a 
d-dimensional (hyper-)cubic lattice as s = 0, . . .2d — 1. 
A sequence of k + 1 steps is then encoded as S = 
(s-k, ■■■So) = (s, so). By P/v, m (S) we denote the sta- 
tistical weight of all iV-step chains in an unbiased sam- 
ple which had followed the sequence S during steps 
TV — m — k, . . . ,N — m. This can be estimated either 
in a previous test run or during the present run. The 
ideal bias in fc-step markovian anticipation is 



2d-l 



p(s |s) = P/v, m (s, S )/ ^ P N,m{s,s' ) 



(8) 



s '=a 



with N 3> to 3> 1 (we used to = 150, and averaged 
fW,m(S) over all N > 300). Thus a step so is made 
more often if this step is anticipated to be more success- 
ful in the far (m steps ahead) future, based on previous 
experience. Eq.(^) should not be used, of course, for 
the very first steps of the chain. There the probabilities 
-Pjv,m(S) are not appropriate. Our remedy is ad hoc but 
efficient: When determining the bias for the n-th step 
in a chain, we replaced fW,m(S) by Pjy, m (S) + const/n, 
with const w 20. The markovian anticipation bias is of 
course compensated by a weight factor oc l/p(so|s), to 
guarantee correct sampling. 

In our actual simulations we used semi-infinite cylin- 
ders. Knowing that the walks anyhow have to grow cither 



parallel or antiparallel to the cylinder axis, we started 
them at height = 1 and put an absorbing barrier at height 
= 0, so that they had to grow into the positive h direc- 
tion. The ratios on the right hand side of eq. (g|) depend 
then on the absolute orientation of the steps and on L. 
Obviously the longitudinal bias is larger for smaller L. 
Most of our simulations were done for the square lattice. 
Results for L = 8 are shown in fig. 1 . There, sequences 
of 9 steps are encoded by integers from to 4 9 — 1, with 
the last step giving the most significant digits. The up- 
ward direction is s = 1, downward is s — 3, and the 
directions perpendicular to the cylinder axis are and 
2. We see that there is both a strong anisotropy, and a 
strong dependence on the shape of the last part of the 
chain. The former means that downward steps are likely 
to be less efficient, while the latter corresponds to the 
fact that a strongly curling walk will have problems to 
be continued. For large L we used isotropic markovian 
anticipation with k = 11. It gave roughly one order of 
magnitude improvement in speed over plain PERM. For 
small L, anisotropic anticipation with k — 8 gave up to 
one additional order of magnitude further improvement. 




FIG. 1. Histogram ratio p(so|s) used for biasing square 
lattice SAW's on cylinders with L — 8, based on sequences 
of length 9. Each sequence S is encoded by a number 
between and 4 9 - 1 = 262143. The probability with 
which the next step is taken is p(so|s), up to corrections 
for the very first steps. On all points, statistical errors 
are < 0.01. The four heavy dots correspond to four se- 
quences S = s, s + 4 8 , s + 2 x 4 s , and s + 3 x 4 s . The 
chain segments corresponding to these sequences are also 
indicated. The probabilities to continue right (r), up («), 
left (I), and down (d), after the segment ururdruu, are 
p(r) : p(u) : p(l) : p(d) = 0.28 : 0.52 : 0.20 : 0. 

We also performed less extensive simulations on the 
triangular lattice, in order to test for universality. There 
we have 6 local directions instead of 4, so that we could 
not use as long memories. We used only isotropic antici- 
pation there, with k = 7. 



2 



In all simulations chain lengths were such that h/L > 
80. Results for the connectivity constant are shown in 
fig. 2 where we plotted \icx> — fi(L) with (j.^ = 2.638159 
p5[ resp. (loo — 4.150795 |17j for the square resp. tri- 
angular lattice. We see perfect agreement with the the- 
oretical prediction indicated by the slope of the straight 
lines. Equation (||) was tested by first making linear fits 
in plots of h versus N for fixed L (to avoid biases from 
the chain ends), and plotting then the slopes against L on 
a log-log plot. The result is shown in fig. 3; again we see 
perfect agreement with the theoretical prediction. Com- 
bining the results from figs. 2 and 3, we find that the 
amplitude B = 0.675 ± 0.002 is indeed universal (the er- 
ror is just a rough but conservative estimate, as are also 
the following error estimates). 
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FIG. 2. Log-log plot of fioo — fi(L) against L. The straight 
lines have slope -4/3, as predicted by eq.(Q), and prefactors 
adjusted to fit the data. Error bars on the data points are 
comparable or smaller than the symbol sizes. 
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FIG. 3. Log-log plot of ho against L, with ho being the 
constant in a fit h = hoN + hi for fixed L. Again, error bars 
on the data points are comparable or smaller than the sym- 
bol sizes. The straight lines have slope —1/3 as predicted by 
eq.M, with prefactors adjusted to fit the data. 



Winding number variances, multiplied by L and di- 
vided by h, are shown in fig. 4. For L < 16 we see fi- 
nite size corrections. Apart from these, all curves co- 
incide within the expected statistical fluctuations, as is 
expected from eq.(||). Again we see that universality is 
satisfied, i.e. the constant A is the same for the square 
and triangular lattices, within the statistical errors. More 
precisely, we find A = 0.475 ± 0.004. 

In order to estimate the universal amplitude D, we 
use the following estimates for the end-to-end distances of 
free chains: lim Ar ^ 00 (i? 2 } /N 2v = 0.7710 ± 0.0004 (16) for 
the square lattice, and 0.711 [|l7j for the triangular lattice. 
Both gave again consistent results, D — 0.9446 ± 0.0006. 

Finally, we consider the number of parallel contacts 
(measured only on the square lattice). A plot analogous 
to fig. 4 shows a strong L dependence, but would not allow 
to estimate this dependence precisely. More instructive 
is a plot analogous to fig. 3. For each L, we first extracted 
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FIG. 4. Plot of L(W 2 )/h against h/L. The horizontal line 
is our best estimate for the constant A. Statistical errors are 
typically of the size of the visible fluctuations. 
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FIG. 5. Log-log plot of ao against L, where m = aoN + ai 
for fixed L. According to the prediction of Ref. one 
should have ai ~ 1/L 1 ' 333 . The straight line is const / L 2 25 . 
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the coefficient ao in a fit nil = ciqN + a%. In fig. 5 we then 
plot ao against L on double-logarithmic scale. From this 
figure we see that nil ~ NL 



-2.25 



nii//i 



-1.92±0.03 



(9) 



Again, the error estimate is just an educated guess. In- 
deed, our data show a slight downward curvature, sug- 
gesting that the true exponent might be closer to —2.0. 
In any case, eq.(6) is clearly ruled out. 

Qualitatively this is indeed not surprising. It means 
that even if the walk winds once around the cylinder, 
the chance that it touches itself is much less than 1 and 
decreases with L. This agrees qualitatively with the be- 
havior of a SAW confined to a strip of width L between 
two repelling walls. In that case the density of monomers 
at a distance z from a wall increases as z x l v fisj] , implying 
that the average distance between two contacts with the 
wall is ~ L 2 . In the present case the previously placed 
part of the SAW acts like a wall, and during one wrap- 
ping around the cylinder the SAW fills a layer of typical 
width L in the longitudinal direction. We thus expect 
in the present case that the number of parallel contacts 
per wrapping is ~ 1/L. This argument is of course not 
rigorous, but it agrees within three standard deviations 
with our numerical result. 

Finally, let us discuss the universal constant A. 
Naively, one could try to estimate it as follows. In planar 
geometry, it is known that the variance of the winding 
number around either end point is fl9lpi 



{W 2 ) = (27r 2 )- 1 lnA 



TV 



(plane) (10) 



(notice that we measure windings in units of 2tt). When 
the plane (punctuated at one of the end points of the 
chain) is mapped onto the cylinder by means of eq.(Q), 
winding around this end point is mapped precisely onto 
wrapping around the cylinder. Taking into account that 
In TV = const + v^ 1 In \z\ —>■ const + 2irh/vL, we would 
predict (W 2 ) — const + (\/irv)h/ L, i.e. A = 1/tw — 
0.4244. This disagrees with our measured value by more 
than 10 standard deviations, and seems definitely ruled 
out. To explain this discrepancy, we notice that confor- 
mal invariance holds only for the canonical (fixed fugacity 
x) ensemble exactly at the critical point (in this ensemble 
TV fluctuates of course, so that we should actually write 
(W 2 J\nN) -► (27T 2 )- 1 for x -> x c = 1/fioo, instead of 
cq.(jl(])). For any finite L, our results for the cylinder 
hold, in contrast, for x c (L) = l/u(L) > x e , since only 
there (TV) — * oo. Therefore, conformal invariance strictly 
spoken doesn't give a prediction for A. 



In summary, we have shown by simulations of very 
long chains on cylinders that universality holds for 2-d 
SAWs, in contrast to recent claims. Our conclusion is 
based on the fact that parallel contacts are much more 
rare than predicted in Q. This suggests that the corre- 
sponding operator, which was predicted to be marginal, 
could play an irrelevant role, so that universality should 
hold for 2-d SAWs. While counting parallel contacts 
with sufficient statistics would have been virtually im- 
possible in planar geometry, it is feasible in cylinder ge- 
ometry. At the same time our results on the constant 
A should serve as a warning that conformal invariance 
should not be applied too naively. 

We thank Gerard Barkema, John Cardy, and Erich 
Eisenriegler for very useful discussions, and the latter 
also for carefully reading the manuscript. 
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